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I. INTRODUCTION 



Since the discovery of high temperature superconductors, much work from both experi- 
mental and theoretical sides has suggested that Fermi-liquid theory does not provide a good 
description of the remarkable normal state properties of these materials. In particular, An- 
derson PJ has proposed that a better framework is to be found in some two dimensional 
generalization of the Luttinger liquid, which has been identified by Haldane 0] as a low en- 
ergy effective theory for many interacting-Fermion systems in one dimension. For example, 
recent measurements of the temperature dependence of the Hall angle in high T c cuprates 
have been interpreted along this line ||. However, it has not yet been possible to derive 
such a picture from a first principles calculation. A suggestive investigation of two-particle 
scattering in two dimensions in the presence of the Fermi sea has lead Anderson to claim 
that a finite phase shift and a singular interaction function f^y exist even in the dilute 
limit . But if the same particle-particle ladder is used in a more conventional many-body 
calculation, it has been shown by Engelbrecht and Randeria |J and Fukuyama et al. || 
that the imaginary part of the self-energy is not singular enough in this T-matrix approxi- 
mation to avoid Fermi-liquid behavior. In fact, the two approaches may not contradict each 
other. Anderson suggests that an extra particle added to the system induces a shift of all 
the momenta of the particles already present in a similar way as the deep level in the X-ray 
edge problem. The present paper is an attempt to connect this physical picture to some 
microscopic perturbative calculations for the repulsive Hubbard model. 

In our approach, we have been guided by the idea already stressed in the work by 
Abrikosov on the Kondo problem 0, and Nozieres et al. on the X-ray edge problem ||, 
namely that both particle and hole scattering channels with the local impurity lead to similar 
singularities which may cancel partially. It is then crucial to treat both in an unbiased 
way. This is for instance a key feature of the parquet summation, which was developed in 
references ]7J and [Q. However, this method is not very easy to implement for many body 
systems, where the angular dependence of interaction vertices is in general not described by 
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a small set of relevant parameters (except in one dimension). A simplification occurs if we 
work in a truncated Hilbert space with only one particle-hole excitation on top of a local 
potential (X-ray edge problem), a spin (Kondo problem), or an added electron (Hubbard 
model). This approach is variational in spirit and amounts to treating, exactly, three- 
particle correlations amongst the particle-hole pair and the analogue of a local scatterer. It 
has been used in the context of Kondo problems to study various properties |]9|,[T0|| , and for 



the Hubbard model in the special case of a single overturned spin [fLlJl . It has even been 
suggested that three-particle correlations between a spin flip and two holes may provide a 
way to reduce the ground state magnetization away from its maximal value in the infinite- 



U Hubbard model at any density (T^|. Three-body correlations are also investigated as a 



possible microscopic mechanism leading to a marginal Fermi-liquid |T3| 



In this paper then, we shall concentrate on the effect of three-particle correlations on 
the self-energy of an electron. We consider a variational state which consists of one electron 
plus at most one extra particle-hole pair excitation on top of an otherwise rigid Fermi sea. 
The scattering in the particle-particle and particle-hole channels are treated on an equal 
footing. In this sense we attempt to go beyond references []5] and || who considered only 
two-particle scattering. We shall be interested in seeing whether the Fermi-liquid starting 
point is valid or not. Exact solutions are available for the Luttinger model and the X-ray 
edge problem and we have decided to compare the electron's Green function in the three- 
body approximation to that of the exact solution for these systems. In section 2 we shall 



discuss the formalism leading to the Faddeev integral equations |14] . We have re-formulated 
these in order to yield directly the self-energy. Sections 3 and 4 will describe results for the 
one and two branch Luttinger models respectively and section 5 will describe results for the 
X-ray problem. 

II. FADDEEV EQUATIONS 

We begin with the definition of the electron Green function, 
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G(q,co) = (N | c^[u-H + i8]- 1 cL + cL[uj + H -i8]- x c^ \ N) (2.1) 



where the Hamiltonian is 



# = E e kcLc kCT + E c k+ P T c kT c k'-p| c k'j. - • (2.2) 

k JVs k,k',p 



We shall consider a variational state where | N) is the non-interacting Fermi sea and the 
Hamiltonian is allowed to create at most one excited particle-hole pair. The corresponding 
diagrams contributing to the Green function are shown in Fig. [I]. We take q > kp and Eq 
is chosen so that (N \ H | N) = 0. In this case, defining 

<P = (N | c qT [uj - H + iS\ I N) (2.3) 

and 

$(k, k') ee (N | c qT [^ - H + z^-^^c^Ck+k^^ I iV> , (2.4) 

the equations of motion of these Green functions truncating terms involving more than three 
particles are 



N 

iVs k,k' 



and 



[Co - e(k, k')]$ (k, k') - ^0 - ^ E $ ( k "- k + k ' - k ") + w E $ ( k "- k ') = (2-6) 

-''s s k" s k" 

where Co = to — U(NJN S ) and e(k, k') = ek + e k ' — e k+ k'- q . The momentum restrictions in 

the sum are such that the arguments of $(k 1; k 2 ) always satisfy kx ^ q; k\^k<i > kp and 
|ki + k 2 — q| < kp. 

We shall find it convenient to define the following functions: 

^(KHiVT'E^K-k') (2.7) 
k' 

and 



4 



j 2 (k) = n; 1 J2HV,k) (2.8) 

k' 

which correspond roughly to T-matrices for the particle-particle and particle-hole scattering 
channels respectively. From the equations of motion for we may see that J\ is related to 
the self-energy by 

S(q,^) = t/ 2 ^J 1 (K) . (2.9) 

K 

From the equation of motion for $ we see that J\ and J 2 satisfy the coupled integral 
equations 

7(K) N-i *UK) U 1 V J 2 (K-kQ 

Jl(K) "^ l — UFi (K) ~ JsT s l — UFjK) \u- e(K K - V) {2A0) 



and 



- L ^(K) U 1 ^ Ji(K + k') 



^-^T^hm+w.TTuFmV^^) ■ (2 - n) 

In the summations the same momentum restrictions as the equations of motion Eqs. (|2.5|) 



and ( |2.6| ) are in effect. F\ and F 2 are defined by 



^(K)^,-'E a _ £(k ,; K _ k0 (2-12) 



and 



f 2 (k) = n; 1 Y,~~ — krJcS ■ ( 2 - 13 ) 

Xf uj- e(k',K) 

The integral equations can be combined into a single one for J\ which is the final equation 
to be solved. 



<A(K) — ]^ 1 _[/^ 1 (k) Sk' (<i- £ (k',K-k'))(l+l/i ;, 2(K-k')) 

U 2 1 V V Ji(k') 

TVf l-UFi(K) ^k" Z^k' (aj_ e (k",K-k"))(l+i7i ; 2(K-k"))(w-e(k'-K+k" ) K-k")) ' 



The technique of truncating the Green function equations of motion at the four point 
level and finding closed equations for them is not limited to well defined particles and is not 
new. It has been applied |L5| to Green functions for the Hubbard 'X' operators [|16[ . 



III. ONE BRANCH LUTTINGER MODEL 



In this section we present the exact analytical solution of the integral equation for a one 
branch Luttinger model (i.e. with right movers only). The total momentum of the three 
particles is q > kp and the energy dispersion is defined by = vp(k-kp). The simplification 
which allows an analytical solution is the fact that e(k, k') = e& + — €k+k'- q = v F (q — kp)- 
Thus we have 



F 1 (K] 



K -2k F 



u> — vp(q — kp) 



and 



F 2 (K) 



q-K 



to — Vp(q — /cp) 



Defining q = q — kp, 



= (u - v F q) l Ji{k) 

fi(k) = (u-v F q)- 1 [l-UF 1 (k)]- 1 
f 2 (k - k") = (u- v F q)- l [l + UF 2 (k - A;")]- 1 



(3.1) 



(3.2) 



(3.3) 
(3.4) 
(3.5) 



we absorb a factor of N s into J\ and go to the continuum obtaining 

fk—kp 



J 1 (k)=f 1 (k) [ ' dk"f 2 (k-k') 
JkF 



rq+kp 

-U 2 dk' 3l (k') 
Jk-k" ■ 



'+kp 



(3.6) 



There is, in fact, no small parameter in this integral equation, even when U is small. 
Thus it cannot be solved perturbatively. This can be demonstrated explicitly by defining 
j(k) = qU 2 ji{k) and x = (u — VFq) /Uq. Rescaling all momenta with respect to q and setting 
kp = we have 

1 r k ..„ 1 



dk"- 



- / dk'j(k') 
Jk-k" 



(3.7) 



x — k Jo x + 1 — k + k" 
Changing the order of integration allows us to perform one of the nested integrals. Using 



k r l 

dk" / dk' 

Jk-k" 



k rk pi rk 

dk' / dk" + dk' dk" 

Jk-k" Jk Jo 



(3.8) 



we find 



x—k 



X 



(3.9) 



The solution proceeds by conversion into a differential equation. Multiplying by x — k 
and differentiating with respect to k we have 



dk'j(k' 



d 1 1 r 1 

^r\(x-k)j(k)] = r r / " y:/,A 

dfc LV m ;J z+l-Jfc x + l-kJk 



Multiplying by x + 1 — and differentiating, 



(3.10) 



[( x + 1 _ _ fc)/(fc) - = j(fc) 



(3.11) 



This can be integrated because the integral equation gives us the two boundary conditions 



j(0) = 



(3.12) 



and 



x{x + l)j'(0) = 1 - / dAtf(fc) = 1 - (Z(x)/xUq) 



(3.13) 



Integration of Eq. (|3.11| ) results in an equation for j(k) in terms of x,k, and j'(0). Integrating 
over k from k = to k = 1 eliminates j'(fc) and and results in an equation for the self-energy 
in terms of x whose solution is (adding in the imaginary parts now) 



= xUqA/(l + A) 



(3.14) 



where 



A = x 2 In 



x 



1-x 2 



1 - Z7rx|x|0(a; + 1)6(1 - x) 



(3.15) 



Substituting back into the electron Green function gives the very simple result 



G(x) = xln 



x 

1-x 2 



- iir\x\6(x + 1)6(1 -x) 



(3.16) 
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where 9(x) is the usual step function. Re- inserting the original units, 



ImG(c<;, q) 



lo — Vpq 



Vpq — Uq < lu < v F q + Uq 



(3.17) 



Uq 



The exact Green function is IT? 



ImG(u), q) = (Uq) 2 - (u - v F qf 



1/2 



Vpq — Uq < uj < vpq + Uq 



(3.18) 



Our result is compared to the exact result in Fig. || The spectral weight from the three- 
body approximation is non-zero over the same energy range as the exact solution and has 
roughly the same shape. Although the two have maxima in the same place, the three-body 
approximation is not able to reproduce the square root singularities of the exact solution. 
Nevertheless spin-charge separation is observed. The charge velocity v c = vp + U and the 
spin velocity v s = Vp — U agree with the exact result. From the integral equations one can see 
that the 'charge' peak comes from [1 — UFi}^ 1 and the 'spin' peak comes from [1 + UF-^" 1 . 

It should be remarked that the one branch model has the same ground state as non- 
interacting Fermions. Thus spin-charge separation can co-exist with a Fermi-liquid. The 
range over which Im G is non-zero is proportional to q — kp and this goes to zero as q — > kp. 
Thus the three-body approximation is self-consistent with the assumption of a rigid Fermi 
surface background. This can be understood by noting that it is impossible to create particle- 
hole excitations from the ground state and conserve momentum. It is this feature of the 
one-branch model which renders the three-body approximation tractable. Below, in the two 
branch model, this phenomenon does not occur. 

Finally we would like to emphasize the non-perturbative nature of the integral equation. 
This feature may help in shedding light on whether or not the two-dimensional Hubbard 
model has a Fermi-liquid ground state. Initial calculations indicate that [1 — UFi]^ 1 and 
[1 + UF 2 ]~ 1 are less singular in two dimensions and almost certainly will not lead to singular 
behaviour on their own. We hope to discuss this in detail in a future publication. 
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IV. TWO BRANCH LUTTINGER MODEL 



We consider a spinless Luttinger model for which only the interaction between opposite 
branches is retained. Again we begin with a particle on the right moving branch but allow 
the creation of a single particle-hole pair on the left moving branch. In this case the integral 
equation is (absorbing a factor of N s into J\), 

1 D 1 

Jl{K) = l-UF^K) kl ,^ +kp (uj-e(k",K-k"))(l + UF 2 (K-k")) 



U 



2 



1 - U F x (K) 



K D D—kp D+K-k' 

E E + E E 

k'=q-k F k"=K+k F k'=K k"=K+k F 



x (4.!) 

(uj — e(k", K — k"))(l + UF 2 (K — k"))(Co — t(k' — K + k", K — k")) V ; 

where the kinetic energy e(k', k — k') = (2k 1 — q — kp)vF (It doesn't depend on k because 
the total momentum must sum to q.), 

F 1 (k)= J2 [Co-eik^k-kT 1 (4-2) 

k'=k+k F 

and 

F 2 (k)= [ D [u-t(k',k)Y l . (4.3) 

Jk' = -k+q 

D is a momentum cutoff. 

We were not able to solve this analytically, mostly because of the intractability of integrals 
of the form J dzln(z) /(z + a), and therefore attempted a numerical solution by discretization. 
In order to control the divergences of the free particle poles we introduced artificial widths. 
That is, we replaced 

1 1 



(4.4) 



uj — q uj — q + iS 

We checked that our results did not depend significantly on the value of 5. The method of 
numerical solution was checked against the analytical solution, Eq. ( |3.16| ) for the one-branch 
model, to confirm its accuracy. 



The exact Green function for this problem is JT? 



lmG(q,w) = - 2D , v ; (a)2 
x 6(u - vq)>y(a, y^)e~ v -y a S l + 6(-Cb - vq)>y(a, -y+)e y - (-y) 



a-1 



(4.5) 



where y± — (Co ± vq)/2D'v, v 



Vp — U 2 and 




(4.6) 



The parameter D' is a momentum cutoff. In a perturbative expansion of the integral equa- 
tions of the three-body approximation the parameter (Co — vpq) /2Dvp appears. In order 
to compare the three-body results to the exact results we need definite values and it seems 
reasonable to take D' = D for this purpose so that the three-body solution and the exact 
solution have the same dimensionless parameter (Co — vpq)/2Dvp. lmG(q, Co) (for the exact 
and three-body solution) as a function of energy with fixed total momentum are superim- 
posed in Fig. [5| Energies are normalized by qvp. We chose a value U/vp = 0.3 which is large 
enough so that the effects of U are apparent but not so large that the integral equations do 
not converge. The cutoff is (D — kp)/q = 3 and is not a sensitive parameter. 

Some differences are immediately evident. Firstly the exact solution has spectral weight 
at Cj < —vq (which is not shown in the figure). The three-body approximation misses this. In 
order to describe this part of the spectrum we would have to include pre-existing particle and 
hole pair excitations in the N-particle ground state whereas the three-body approximation 
assumes that the background is a rigid Fermi sea. A background of pre-excited particle- 
hole pairs would require consideration of a minimum of five bodies. The second difference 
is that the spectral weight of the exact solution is extremely concentrated just above the 
maximum of the spectral function. That is because, for small U, Im G diverges with a 
negative exponent only slightly greater than —1. The exact solution displays an interesting 
characteristic. For energies \oj\ < vq the spectral weight is exactly zero. We do not know 
how to explain this in terms of the underlying electrons. This feature is not present in the 
three-body approximation. Nevertheless the spectral function in that case is asymmetric. 
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In Fig. ^| the reflection of the spectral function about its maximum is plotted to bring out 
the asymmetry. We see that the asymmetry is in the right direction. That is, the spectral 
weight is higher for energies above the maximum. By looking at the numerical solution 
integral equations it is possible to tell that a pole in [1 + F 2 ]~ 1 (the particle - hole scattering 
channel) is the cause of this asymmetry. 

Another feature which is reproduced qualitatively is the negative energy shift of the 
peak position (real part of the self-energy). In the exact solution this manifests itself in a 
downward renormalization of the Fermi velocity v = \Jv 2 F — U 2 . In Fig. [| the shifts in the 
maxima of Im G as a function of U are plotted. It is interesting to note that the three- 
body approximation follows the square root behaviour of the exact solution quite accurately 
except for a factor of two. 

To summarize, the three-body approach reproduces qualitatively the real part of the 
self-energy and the asymmetry of the spectral function. Although it is not shown in Fig. |3|, 
in the three-body approximation the spectral weight of an electron with momentum q > kp 
stretches into negative energy, with a width proportional to q — kp. However, the integrated 
spectral weight for negative energies seems to be independent of q — kp at small q — kp and 
amounts to a few percent of the total spectral weight (for the given choice of parameters): 
comparable to that of the exact solution. This flow of spectral weight can be interpreted 
by saying that the three-body approximation is inducing some correlations which tend to 
improve the trial ground state. This confirms the idea that correcting the N-particle ground 
state by coherent particle-hole pair excitations would yield a Green function more like the 
exact solution. 

V. X-RAY PROBLEM 

As discussed in the introduction, our goal was to set up an approximation scheme which 
treats particle-particle and particle-hole interactions on the same footing. Formally, this idea 
was first discussed in condensed matter physics for the Kondo problem |7| , and for the simpler 
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X-ray edge problem ||. In this section, we present the results of our simple three-body type 
of approach for the X-ray edge problem. At this point, it would be helpful to characterize 
better the expected limitations of our method. The quantity of interest here is the overlap 
between the wave function of the conduction electron system, after the sudden turning on of 
a localized potential, and the unperturbed Fermi sea. Going beyond the first calculation by 
Mahan ||18| , Nozieres and coworkers have shown that the deep level Green's function defined 



as above decays as a power- law function of time |19| , with an exponent proportional to the 



square of the phase shift at the Fermi energy due to the local potential. This result has 
been interpreted in terms of Anderson's orthogonality catastrophe pO| , pT| , and rederived in 
the simpler language of Tomonaga bosons [ 22| . As stressed by many authors, the infrared 



singularity arises from the excitation of a logarithmically divergent number of Tomonaga 
bosons. However, the average number of such emitted bosons increases only logarithmically 
with time, which might enable us to restrict ourselves to the subspace with no more than one 
excited boson, at least for times shorter than -j^exp (W/2V 2 ), where W is the conduction 
electron bandwidth and V is the strength of the localized potential. 
The Hamiltonian for the X-ray problem is 

H = £ bb j + J2 ek4 a k - E + bb ] E ^k',k4' a k • (5.1) 

k kk' 

at is a spinless Fermion annihiliation operator and Sq is the energy of the deep hole. As 



in the Hamiltonian |2j|, Eq is chosen so that the unperturbed Fermi sea with the deep 
level occupied is defined to have zero energy. Let \N) be the unperturbed Fermi sea and 
|0) = b\N) be the unperturbed sea with the deep hole. We shall be interested in the Green 
function 

Q(t > 0) = -t(N\e iH Ve- iHt b\N) = — ze~ i£:o * <0 1 e"^* 1 0> (5.2) 

We again set up the three-body approximation using the equation of motion formalism. 
Defining \t) = exp (— iHt)\0), we make the approximation 

\t)=m0)+ E E $(k,k\t)4a k ,|0> (5.3) 

k>kp k'<kp 
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with the truncated equations of motion 



i-£(t) =S <J)(t) + J2Vk>,kHk,k\t) (5.4) 

k.k' 



and 



<9$ 

i— (k, k\ t) = (S + e k - eu)$(k, k', t) + E fc » V k}k »®{k", k', t) 

- Ek" V k ,, k ^(k, k", t) + VJfe, fc ^(t) . (5.5) 

We proceed by expressing $(k, k', t) in the basis of single-particle eigenstates in the presence 
of scattering V k>k >, 

=£s«^(f)#(A0$(AO ( 5 - 6 ) 

in which case the second equation of motion may be written 

^CW^fc')^(t) = E(^ + 4 + WWW^W+^W • ( 5 - 7 ) 

a,/3 a,/3 

Using the orthogonality of this basis we have 

$a,p(t) = -i fdt'exp \-i{El + E h p + S )(t - 0]^ECW^(^)^ ( 5 - 8 ) 



k,k' 



Combining this with the first equation of motion written in the same basis and defining 
4>(t) = exp (-i£ t)4>(t), 



0(t) =1+1* dt' f dt" £ £ Vk,lVG^(t' - - t"W') (5.9) 



where 



G£J(*) = E e-* ( V/WCl?) (5-10) 



are the particle and hole Green functions in the presence of the deep hole. Upon Fourier 
transforming Eq. (|5.9| ) and solving we find that <f>(uj) — — S(a?)) where the self-energy 
£ is given by 

sH = iEE» / ^^(^^(w-wO (5.ii) 
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We take the continuum limit and treat the deep hole as a point scatterer thus dropping 
momentum dependence in V. The self-energy becomes 

S = 1V 2 ( ^G p (uj')G h (uj - uj') (5.12) 

J 2,71 

where G p ' h are now the on-site Green functions given by 

(l/2Dv F )\n [(uj - v F k F + i5)/(u - v F D + i6)] 
K ] 1 - (V/2Dv F )\n [(uj - v F k F + i6)/(u - v F D + iS)] 1 J 

and 

h (l/2Dv F )\n [(u + v F k F + i5)/(u - v F D + i5)] 

{) l + (V/2Dv F )\n[(u + v F k F + iS)/(uj-v F D + i8)} ' { V 

D is the usual momentum cutoff. 

In Fig. [5] we have plotted the imaginary part of the particle and hole Green functions. 
Note the weak logarithmic singularities which result from single particle bound states. In 
Fig. ^ is plotted the resulting real and imaginary parts of the self-energy. The critical 
features in the latter plot are the sharpness of the cusp in Im£(u;), and the steepness of 
the inflection point of ReS(o;) near uj ~ 0. As this numerical calculation shows, the weak 
logarithmic divergences of the single particle spectral density do not have a great effect on 
the self-energy. 

In view of this, we proceed with an analytical calculation of the quasiparticle residue 
which replaces the spectral functions plotted in Fig. |5]by the constant —(2Dv F )~ 1 multiplied 
by appropriate step functions. For small uj, 

Re£H (V/2Dv F ) 2 [v F (D - k F )ln (^) + v F (D + k F )\n (^) 

^( , F ,'S Wf ) )] ■ ( 5 - 15 ) 

The pole in the deep hole Green function occurs at u p = Re£(u;p) which is 

V 2 



2Dv F 

in the k F « D limit. The residue of the deep level pole is 
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(ln2 + 0{k F /D) 2 ) (5.16) 



= 1 - {V/2Dv F f\n . (5.17) 

These expressions have been quantitatively checked with a numerical calculation. 

This calculation shows that the simple three-body approach does not reproduce the 
power-law singularity of the deep level Green's function at low energy. Instead of a power- 
law decay at long times, we obtain a long-lived level, with a residue close to unity at small 
coupling. This result is by itself not too surprising. It simply confirms the intuitive idea that 
the possibility of exciting many particle-hole pairs should be retained, in order to describe 
the power-law singularities. In fact, this feature is present in the parquet diagram approach 
used by Abrikosov and Nozieres and coworkers |7||§. The similarity between a parquet 
calculation and our simple three-body approach arises because in the former the three- 
particle vertex function involves only successive two-particle interactions. However, since 
the same parquet diagram corresponds to several time orderings for the interaction events, 
the actual number of particle-hole pairs is not fixed to unity. Unfortunately, our simplified 
approach does not seem to provide a way to bypass the cumbersome parquet algebra. The 
main complication encountered upon going from the Faddeev approach, with fixed particle 
number, to the parquet algebra in the full Hilbert space, is reflected by the need to keep 
frequencies as additional integration variables in the internal lines of the graphs. For the X- 
ray edge problem with a separable potential, momentum integrations are staightforward and 
we recover a one-dimensional problem (in the time direction). For more realistic systems, 
such as the two-dimensional Hubbard model, such a simplification is not so obvious. This 
is why our simpler method may still be useful there, since it may for instance be able to 
indicate the presence of spin and charge decoupling. We hope to be able to address this 
issue in a forthcoming investigation. 
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FIGURES 

FIG. 1. Diagrams in the one-electron self-energy accounted for in the Faddeev Equations. As 
opposed to the usual convention the vertices are time ordered from left to right. This prevents a 
violation of the variational assumption of the presence of at most three excitations at one time. 
The middle line is the spin up electron and the outer lines are the spin down particle-hole pair. 

FIG. 2. One branch Luttinger model: comparison of Im G for the exact solution (dashed line) 
and the exact solution of the three-body approximation (solid line) as a function of normalized 
energy x = (Co — vpq)/Uq for fixed total momentum q. 

FIG. 3. Two branch Luttinger model: Im G as a funtion of normalized energy CbjvFq. Solid 
line: three-body approximation; Short dashed line: exact solution; Long dashed line: three-body 
approximation reflected about its maximum. The non-zero part of Im G at negative energies is 
not shown. 

FIG. 4. Two branch Luttinger model: Real part of the self-energy defined as the energy shift 
(relative to and normalized by the non- interacting quasi-particle energy vpq) of the peak of Im G 
as a function of the interaction strength U jvp. Dashed line: exact solution; Diamonds: three-body 
approximation . 

FIG. 5. Spectral density of the on-site particle (solid line) and hole (dashed line) Green 
functions in the presence of the deep hole. The values Hf/D = 0.1 and V/Dvf = 0.2 have been 
used. 

FIG. 6. Self energy of the deep level as a function of frequency using the same parameters as 
in Fig. |5[ Real part: solid line, Imaginary part: dashed line. 
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